Topological Charge Order in a Frustrated XY Model and Related systems 
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We prove the existence of a finite temperature Z2 phase transition for the topological charge order- 
ing within the Fully Frustrated XY Model. Our method enables a proof of the topological charge 
confinement within the conventional XY models from a rather general vista. One of the compli- 
cations that we face is the non-exact equivalence of the continuous (angular) XY model and its 
discrete topological charge dual. In reality, the energy spectra of the various topological sectors are 
highly nested much unlike that suggested by the discrete dual models. We surmount these difficul- 
ties by exploiting the Reflection Positivity symmetry that this periodic flux phase model possesses. 
The techniques introuduced here may also be applied to examine transitions associated with other 
topological defects, e.g. the confinement of disclinations in the isotropic to nematic transition. 



I. INTRODUCTION 

Here we examine the topological charge ordering of XY 
models in the presence of various external magnetic fields 
and prove, by Reflection Positivity, that these charges 
order at finite temperature in the "Fully Frustrated XY 
Model" (wherein half a fluxon threading each plaquette). 

The outline of the paper is as follows: In section (II) 
we define the two dimensional frustrated XY model. In 
section (III), a simple set of general trigonometric dif- 
ference equations and approximate solutions for ground 
state configurations will be given. Much of this section 
is not new but is merely our perspective on the matter. 
We illustrate how in the low- energy sector the continu- 
ous XY model and its discrete topological charge dual are 
almost one and the same. In section(IV), several simple 
exact ground state solutions are given. We then compute 
energy gaps between different topological charge config- 
urations which will serve as a springboard in our proof 
of finite temperature phase transitions. 

Henceforth, we report new results. In section(V) we in- 
vestigate the Fully Frustrated XY model. Employing sec- 
tion (IV), in subsection(V A) we establish the existence 
of energy gaps between the different topological charge 
sectors. We then show, in subsection (VB), that the sys- 
tem is Reflection Positive and consequently prove a fi- 
nite temperature topological charge ordering transition. 
In section(VI) we provide an alternative proof to the or- 
dering transition. We demonstrate, in section(VII), that 
by a minor variation, our proof immediately yields the 
usual finite temperature annihilation of vortex-antivortex 
pairs occuring in the Kosterlitz-Thouless transition. In 
section (VIII) we discuss systems with kaleidoscope pat- 
tern frustrations. In section (IX), we discuss the exten- 
sion of the chiral Z2 symmetry of the Fully Frustrated 
XY Model to a Zk symmetry present for other rational 
frustration / which tends to an 0(2) symmetry for irra- 
tional /. The Z2 symmetry which is broken in the Fully 
Frustrated XY Model might be a hidden degenerate form 



of an 0(2) symmetry. Numerically, both symmetries 
are broken near unison within the Fully Frustrated XY 
Model (at T ~ 1.286J (with J the exchange constant)). 
We suggest an analogy to the magnetic groups appearing 
in the Quantum Hall problem. We end, in section(X), in 
an application of these ideas to other topological defects- 
e.g. disclinations in liquid crystals. 



II. DEFINITION OF THE MODEL 

We consider a classical XY model on an L x L square 
lattice with periodic or open boundary conditions (with 
2L^ or 2L^-|-2L bonds respectively). The two component 
spins {Si} at each lattice site i are labeled by their ori- 
entation {9i} relative to a chosen axis. The Hamiltonian 
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where i,j are nearest neighbor sites on the lattice. Here, 
as well as throughout the paper, we set both the ex- 
change constant J, and the lattice constant to one. Two 
symmetries of this Hamiltonian are evident: 

• A global 0(2) symmetry. The transformation 0i -^ 
Oi + 4> for all sites i with a uniform shift (f) leaves 
the Hamiltonian H trivially invariant. By the Mermin- 
Wagner-Coleman [1,2] theorem, continuous symmetries 
are never broken in two dimensions: at any finite tem- 
perature the local order parameter {6i) =0 for all twice 
differentiable, rotationally symmetric, momentum space 
interaction kernels [3]. In this two dimensional system, 
the superfluid phase stiffness exhibits a first order jump 
at a Kosterlitz-Thouless temperature Tkt at which ro- 
tational symmetry is (marginally) "broken" [4], [5]. 

• For the special cases in which the flux (f>/j 
is the same for all plaquettes R, translational invariance 
is also a symmetry of the Hamiltonian. We shall pri- 
marily focus (for the Fully Frustrated XY Model) on 
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the breaking of translational invariance under a shift by 
one lattice constant, 9i -^ 9i+ei for all i, as seen by the 
arrangement of vortices. A critical temperature T^ortex 
may be associated with the breaking of this symmetry. 
Numerically, as the temperature is lowered within the 
FFXY, it is found that \T^ortex-TKT\ < Tkt- Whether 
the two temperatures might be exactly one and the same 
has been a puzzle surrounded by some controversy, albeit 
very promising recent work [6] . 

Eq. (1) describes an array of Joscphson junctions im- 
mersed in an external magnetic field. It may model an 
extreme type II superconductor built of discrete super- 
conducting elements ( "grains" ) at the various lattice sites 
where it is assumed that the modulus of the supercon- 
ducting wave-function on each grain is pinned and that 
only its phase which may vary [7] . The evolution of the 
phases is easily gleaned within a path integral formula- 
tion: an external electromagnetic vector potential A^ (x) 
gives rise to the differential action dSem ~ —^A^dx^^ 
for a displacement dx^ , with e* the effective charge 
and c the speed of light. This leads to a phase shift 
exp[i5'em/^] ^ cxp[iAij] along the link {ij) where 
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It follows that in going around a plaquette the electro- 
magnetic field incurs an additional phase 
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the elegant Aharonov-Bohm effect [8]. The above di- 
rected loop sum may be replaced (by a simple applica- 
tion of Stokes theorem) by the magnetic flux threading 
the plaquette to conform with the conventional experi- 
mental prognosis of the effect. The elementary plaquette 
R is said to be frustrated if the Aharonov-Bohm phase 
associated with it is not congruent to zero. In such an 
instance, the net magnetic flux piercing R is not an inte- 
ger multiple of the elementary fluxon $o- If the directed 
sum is defined as 
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then the "frustration" /^ is not an integer in such a case. 
For the type II superconductor, / = ^^ with a the 
lattice constant, B the magnetic field strength, and $0 ~ 
■|| the elementary magnetic fiuxon (e* — — 2e for the 
Cooper pairs) . All of the above may be formulated in a 
gauge independent manner. We define (pij = 9i — 9j — 
Aij (mod 2tt) such that — tt < (f>ij < tt. The Hamiltonian 
is now a function of {(t>ij}. Unimportant gauge degrees 
freedom have been absorbed. In a system with periodic 
boundary conditions, we are left with N — L"^ constraints 
of net circulation, one for each plaquette R: 
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{ E 't^ij = 2iT{mR - Jr) } 



(5) 



with integer {ruR} [9]. The content of the "vorticity" 
•niR is clear. Let us slowly trace its origin and justify 
the latter name. The vorticity of a plaquette R: mn — 
J2{ij)&R^ij where (pij = 9i - 9j - Aij + 27rmy. The 
vorticity iriR represents the multiple of (27r ) that needs to 
be added to the bare gauge invariant angular differences 
{9i — 9j — Aij) such that — tt < (t)ij < tt. As such, it 
counts the number of times that the XY rotor circulates 
as we go counterclockwise round the plaquette R. We can 
immediately make two statements regarding the effect of 
a spatially uniform frustration /: 

(i) The ground state energies are trivially continuous 
in /. The proof is very simple: By Eq.(l), we note 
that for each given angular configuration {9i}, the en- 
ergy E{{Aij};{9i}) is a continuous function of {^y}; 
each individual gauge link variable Aij appears as an 
argument of single continuous cosine function. Thus 
the ground state energy (i.e. the minimum amongst 
all of these continuous energy curves, each correspond- 
ing to a different angular configuration {9i}), Eground ~ 
infe{i?({Ay}; {^i})} is also a continuous function of 
{Ay} or, in a gauge invariant formulation, the ground 
state energy per link is a continuous function of the frus- 
tration / [10]. This continuity argument is shorter than 
that spanned by an entire earlier Letter [11]. 

(ii) The system is invariant under the transformations 
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for all integers n. Once again, the proof is trivial: 
the inversion of the external magnetic field / -^ — / 
along with 9i -^ —9i leaves the partition function Z ~ 
Hi / <^^i cxp[— /3i/] identically the same as the cosine in- 
teraction term is an even function. The periodicity of the 
cosine secures an additional invariance: f ^ n + f [12]. 
Compounding the last two operations together proves 
the symmetries of Eq.(6). These imply that, as a func- 
tion of external magnetic flux, the ground state energy 
Egroundif) has local extrema at the symmetry points 
/ = 0,±1/2,±1,... . The points / = (mod 1) are 
global minima of Egroundif)- all bonds are saturated at 
their minimal value (—1) within the ground state. 

The symmetries of Eq.(6) have led the / = 1/2 vari- 
ant to be regarded as the "Fully Frustrated" one. After 
all, any larger value of / may be folded back to another 
value of / lying in the interval [0, 1/2]; the model with 
/ = 1/2 is as frustrated as it can possibly get. How- 
ever, in some sense, the "Fully Frustrated XY model" 
is the least frustrated of all. The time reversal symme- 
try (or, more formally, a Reflection Positivity invariance) 
that this model possesses will allow us to prove that it 
displays a Z2 phase transition. 



III. APPROXIMATE GROUND STATES AND 

THE LOW ENERGY DUAL (TOPOLOGICAL 

CHARGE) MODEL 

The current section, unlike others to foUow, is not rig- 
orous. Much of its contents are not new (albeit often de- 
rived by new routes). Its intent is to provide the reader 
with better intuitive grip on the physics of the frustrated 
XY model. Experts may choose to directly study Sec- 
tions IV or V. Enforcing the constraints of given vortic- 
ities, Eqs.(5), the Hamiltonian to be minimized reads 
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where {Aq} are Lagrange multipliers. Unless some of the 
angles lie on the boundaries {4> = ±7r), at the cxtrcma 
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where the two adjacent plaquettes Ra,Ra' share the 
bond (ij). The Lagrange multipliersAyga and AJ g^ re- 
fer to the plaquettes a and a' which share the same bond 
(ij) (one having it appear in a counterclockwise circula- 
tion of a and the other having it in a clockwise reading 
of a'). Examining Fig.(l) we sec that 



sin 001 = A/ 
sin (f>03 = Xiii 



A//, sin(/)04 = 
Xiv, sin(/)02 



A// — A/77 
= A/y — A/. 



(9) 
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FIG. 1. A lattice fragment centered about the origin. 



Summing it all up, 
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(10) 



the sum of all Joscphson currents entering any node 
must vanish. We also arrive at Eq.(lO) when minimizing 
Eq.(l) within a given gauge to find the global ground 
state. Eqs.(7,10) show that current conservation holds 
at extrema of each of the individual topological sectors 
(toi, m2..., niN) (unless the minima lie on the boundaries 
of the topological sectors: i.e. (pij = ±tt). 



If / = and mji = in the region outside a given 
set of vortices then we may work in a gauge in which 



6i — 6j = 6ij. In this gauge, Eq.(lO) reads 
V^(sin[V^0]) - 0, 



(11) 



with V^ denoting the discrete lattice difference along the 
fi direction. This is equivalent to the difference equations 



W9 = sm-^ ([W xA{x)]'' 



The function A{x) must satisfy 
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for all contours T. Here /iy are chosen to lie along each of 
the individual bond directions (ij) in the contour prod- 
uct. In the continuum limit Eq.(ll) becomes 



cos^ -df.d'^eix) = D^e{x) = 0. 



(14) 



Effectively, the phase 9{x) isotropically dilates the met- 
ric. Apart from the special case 9 = ±7r, we obtain the 
Laplace equation. Thus, in the vortex free {m ~ 0), frus- 
tration free (/ = 0), region 9{x) is a harmonic function. 
Let us now return to the discrete lattice equation Eq.(ll) 
and examine its behavior in the linear regime far away 
from the vortices when the angular differences are small 
we can replace Eq.(ll) by V^V^0 = 0. Viewing the 
Joscphson current as a fictitious magnetic field B = W9 
(with V the discrete lattice gradient), the condition for 
a minimum becomes V • ,8 = at all lattice sites. On 
the other hand, the boundary conditions (recall that we 
focused on a region in which mjj — fn = such that a 
decomposition (j)ij — 9i — 9jis possible) imply that 
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where m„et and /„et are respectively the net topologi- 
cal charge and external magnetic flux/(27r) in the region 
bounded by the contour C. Choosing a gauge for the 
pscudo magnetic field B = ^7 x A such that V • ^ = we 
obtain an equation similar, in the thermodynamic limit, 
to the lattice version of the standard Biot-Savart law 



Aix) = 27rJ2[ms'-fs']Gix,: 



G{x, x) 



exp[zfc • (x — x')] 



B.Z. 



d^k 



(27r)2 4- 2cosfci - 2cosA:2' 



(16) 



where m^ and fg are the local vorticity and frustration in 
the plaquette x. The sum over x' spans the entire lattice. 
(In higher dimensions the denominator readily general- 
izes to 2 ^j^^(l — cos A:/).) Employing the "Biot-Savart" 



law and the magnetic analogy, boundary effects may be 
included by extending the integral to include images of 
the singularities. The ground state energy given a certain 
distribution of {mg} and {fs} outside our domain is 



Aground — ^ / ^ y^x Jx\^\^ ^ -^ )\P^x' Jx'\ 



(17) 



In the small k (continuum) limit G{x — x') becomes 
logarithmic. The resulting expression becomes identical, 
for all practical purposes, to the Villain model. In this 
limit, the energy of interaction between two vortices is 



E = 2'Kqiq2\n{a/ri2). 



(18) 



where qi = (m^ — fj). For the time being, wc have in- 
serted a length "a" in order to make the argument of the 
logarithm dimensionless. When supercurrent screening is 
introduced, "a" naturally becomes the screening length 
in the short distance limit. Thus, in this unscreened case, 
"a" may be regarded as the system size: a = 0{L) [13]. 
In the continuum limit lines connecting opposite sign 
vortices become cuts in the complex plane with branch 
points at the positions of the vertices. A solution to 
the Laplace equation (Eq.(14)) jumping by 2Trqi across 
the branch cut while traversing each pole (vortex) Zi — 
Xi 1 +ixi 2 (where {xi i,Xi 2) are the vortex coordinates 
within the plane) is the well known [14] 



d{z) = ^ qi/m{ln(z - Zi)}. 



(19) 



In the previous (magnetic) analogy the singular branch 
cut takes on the role of a Dirac like string. If we have 
a pair of equal and opposite vortices qi = —q2, then 
9/qi simply becomes the "angle of site" of the segment 
connecting zi and Z2 as viewed from the point z. 9 may 
be written as the phase of wave-function (which may be 
superficially viewed as a Laughlin-like wavefunction as a 
function of only one of the "electronic" coordinates) with 
the odd integer power replaced by q [15]: 



^l;{z)=X{{z-z,f 



(20) 



The phase of this wavefunction 9 — Im{\nil){z)}. 
The continuum limit unfrustrated XY Hamiltonian, 



H 



(fx{vef 



hi 



-- - / (fxVIm{\ml;{z)} ■ VIm{lntp{z)} 
(fxVRe{\n{z - Zi)} ■ Vi?e{ln(z - Zj)}, 



where we employ the Cauchy-Riemann relations 

daU = €a(3d(lV =^ (Vu)^ = {^vf (21) 



(with £12 = —£21 — 1) for the complex function F{z) = 
\'n'ip{z) = u + iv (where the squared gradients are l/r^). 
Integrating by parts, 

H = --^qiqj (fxln\z - Zi\\/^ln\z - Zj\. (22) 

By the 27r discontinuity around a logarithmic pole, 

V^ In \z — Zj\ ~ 2'!rS{z — Zj). Thus, 

H = -TT^ qiq-j In |zi - Zj \ 
= TT 2^ <liqj lii(a/^ij) + const (23) 



with an arbitrary a and inter-vortex separation r^ = 
\zi — Zj\. This swift independent derivation is merely a 
slight twist on the treatment of Itzykson and Drouffe [17]. 
Note that the final result is identical to Eq.(18) derived 
by a seemingly alternate route. Such a plasma like inter- 
action is also present in the Boltzmann weights of Quan- 
tum Hall Laughlin states. Along this route, the standard 
Villain model duality [18] acquires a new meaning- an 
interchange between the real and imaginary parts of a 
complex function F(z). Here the nature of the duality 
is manifest: the imaginary part of F(z) leads to a de- 
scription in terms of the angular configuration of the XY 
rotors. The real part of F{z) leads to a description in 
terms of the topological charges of the defects (when the 
frustration / = 0). Formally, we linked the two dual pic- 
tures by employing the Cauchy-Riemann relations. We 
note that for / = 1/2, the Z2 transition temperature 
T'uortex corresponds to the ordering of charges in the dual 
model, whereas Tkt corresponds to the critical temper- 
ature in the original angular variables. Number (charge) 
and phase variables appear in many instances (e.g. in 
the triangular XXZ models displaying both Z2 and XY 
transitions with nearly equivalent temperatures [19]). 

We now provide a new geometrical interpretation of 
the frustrations. Forgetting about boundary conditions, 
let us arrange the lattice sites on a sphere of radius R (or 
many fractured portions of such spheres pasted together) . 
The sum of the angles of a spherical quadrilateral = 2tt 
+ Area/i?^. It follows that in this naive mapping, we can 
identify a uniform frustration / = (Area of the plaquette) 
/ R^. The frustrated XY Hamiltonian is now 



H 



{ij) 



o S-i 



(24) 



where o denotes the scalar product between the vector 
resulting from parallel transporting the spin Si to site j 
and the spin Sj already situated at site j. This Hamil- 
tonian is translationally invariant when / = const and 
may be diagonalized in Fourier space [22] . By this crude 
analysis, we would expect the ground state modulation 



length (for a constant value of /) to roughly scale as the 
radius R ~ f~^^^ for small fractions / such that one 
fluxon pierces the surface of the sphere (giving rise to a 
trivial Aharonov Bohni phase). This conforms with nu- 
merical observations for f — l/q with an integer q ^ 1: 
each fluxon (or vortex) occupies an area of size 1// [20]. 
More conventionally, requiring overall charge neutral- 
ity in the Villain model (or for our pseudo-currents {qa} 
that interact Coulombically) yields the same result. By 
net charge neutrality, the background charge [— ^^ /„] 
should cancel against [^^ ma] originating from the vor- 
ticities. Once again, from this "angle" as well, the den- 
sity of vortices should scale as /^/^ for a small uniform 
frustration /. In any system having periodic boundary 
conditions, the net topological charge must vanish. It 
follows that any system having ferromagnetic boundary 
conditions cannot enclose a net nonzero vorticity. All 
vortices must appear in vortex-antivortcx pairs. On a 
less precise level, within the approximate dual model, 
this follows from a divergent Coulomb penalty for any 
lone unpaired vortex. Let us define a radius r vortex in 
the unfrustrated problem as one in which 



(y>er 



= 27r?7i 



R, 



(25) 



with rriji = 1 is satisfied for all contours F which encircle 
the plaquette i? up to a maximal distance r away from 
it. Employing the Euler-Maclaurin summation formula, 
we find that for large r, the energy of a radius r vortex 
is bounded from below by that of a vortex in the con- 
tinuum limit where an integration may be performed. In 
the large r limit the energy may bounded from below 
by Ki ln(r/ri) where, in this unfrustrated case, the con- 
stants < Ki < IT and ri > 1. "Phase space" volume 
calculations may also be enacted for these systems [23]. 



IV. EXACT FLUX PATTERNS 

All that follows in this section and hereafter is exact 
and is not based on any of the approximations in sec- 
tion (III). The periodicity (if any) of the ground state 
phase pattern {9i\ depends upon the gauge. Even in the 
unfrustrated (/ = 0) case a gauge with irrational {^y} 
may be chosen such that the phases are never periodic. 
As throughout most of this paper we will attempt to focus 
on the gauge invariant bond angles {0ij}- A periodicity 
in {0ij} implies a periodicity in the physically measur- 
able Josephson currents (sin(/)y). It is easy to determine, 
analytically, the ground state of small fragments of the 
lattice. As before, we enforce the constraints via La- 
grange multipliers; now we leave the extrema equations 
in their plain form and resort to solving a few trigonomet- 
ric equations and, in the aftermath, compare contenders 
amongst the resulting constrained extrema (along with 



various states that lie on the boundaries of {li-bonds ® (f'ij } 
of the different topological charge sectors: scenarios in 
which a few of the bond angles are <j>ij ~ ±7r). 

Let us start by considering a single plaquette with 
open boundary conditions. If this plaquette is unfrus- 
trated, the ground state is trivially a ferromagnet: all 
(j>ij — (this state lies in the m = sector). If this 
plaquette is threaded by half a fluxon (/ = 1/2) then 
within the ground state the gauge invariant angles are 
all the same 4>i2 = ^23 = ^34 = 4>ii = ±7r/4 where 
i ~ 1,2, 3, and 4 are the vertices of the plaquette. (Note 
that the directions in these bond angles are important 
(j)ij = —(f>ji). If all angles are 7r/4 then the plaquette 
has topological charge q = 1/2 or a vorticity to = 1; if 
all bond angles arc — 7r/4 then the plaquette charge is 
q = —1/2 and the lies in the m = sector. Other topo- 
logical charge sectors are elevated relative to the ground 
state by a finite gap. For instance, the minimum energy 
for the unfrustrated plaquette (/ ~ 0) subject to the 
constraint to = 1 is e™~^ = —1/2 (with a configuration 
having (j)ij — 7r/3 for three bonds and with the remaining 
bond on the boundary of the angular zone {(jjij = tt)) lies 
above the ground state e™~f = —4. Thus, each state in 
the m ~ I sector is elevated by at least 7/2 (J- if units 
are restored) relative to the ground state [16]. Similar 
bounds may be generated in frustrated plaquettes. These 
finite energy gaps will allow us, later on, to argue for fi- 
nite temperature vortex ordering for certain values of / 
and to prove the existence of a Z2 phase transition for 
the Fully Frustrated XY model observed numerically in 
[24]. Note that here the ground states computed with 
open boundary conditions can trivially tile the plane. 
One of the complications of the continuous XY spectrum 
is that the various topological charge configurations are 
nested. The energies of the single unfrustrated plaque- 
tte span the region —4 < e(m = 0) < 4 when m — 
and —1/2 < e(|TO.| = 1) < 4 when a vortex/antivortex 
is present. The energy of a vortex free plaquette may 
be higher than that of a vortex-full plaquette! The nest- 
ing of the energy spectrum just found hints that the XY 
model cannot be probed by merely replacing its plaque- 
ttes by topological charges. This subtle non equivalence 
between the continuous model and its discrete counter- 
part is not heavily emphasized in the literature. The dual 
model (derived by various approximations for the low en- 
ergy sector in Eq.(17) and Eq.(23)) is a useful tool but it 
is merely an approximation albeit an amazingly beauti- 
ful and powerful one. The complication arising from the 
non-equivalence of these models will be one of major ob- 
stacles that we will need to surmount. As noted, the only 
ground states for a uniform / = 1/2 are the two checker- 
board patterns where <t>ij£R = dnr/A for all plaquettes R. 
The proof is quite simple: the energy 



-S^^^E'^^-oE ™in{e„} 



(26) 



where e^ is the energy of each individual plaqucttc. Thus 
if a configuration in which for each plaquette ea — emin 
exists then that configuration is a minimum energy con- 
figuration. For each individual plaquette i?Q., if Cq = Cmin 



then all 



HjGa 



±7r/4. Now, a bond {ij) is common to 



two adjacent plaquettes. If its value, as read counter- 
clockwise from the center of one plaquette, is 7r/4 then 
its value as read counterclockwise from the center of the 
adjacent plaquette is — 7r/4. If we start from one plaque- 
tte in which all bond angles are 7r/4 and then add an 
adjacent plaquette and require that the in the new ad- 
ditional plaquette all bond angles will be of equal value 
which is set to either 7r/4 or — 7r/4 then we will triv- 
ially see the bond common to both plaquettes will force 
all bonds in the new plaquette to be equal to — 7r/4. 
We can then add a new additional plaquette to the last 
one and repeat the argument. Thus if within the "seed 
plaquette" (the plaquette situated at the origin of the 
dual lattice) all gauge invariant bond angles (l){ij) = 7r/4 
then any other plaquette lying at a point x = (a;i,a;2) 
will also saturate the minimum energy bound with all 
bond angles set to [(— l)('^i+^^^]7r/4. We could have sim- 
ilarly started with a plaquette at the origin in which all 
bond angles were —n/A. These are the only two possible 
"seeds" and the resultant configurations are the only two 
global ground states. Within these two ground states 
we find a checkerboard pattern of topological charges: 
If, in a given plaquette R, all the gauge invariant bond 
angles (/)ij = 7r/4 then the vorticity mjj — 1; when all 
(l)ij — — 7r/4 there is no vorticity {mji = 0). The one and 
zero values of the local vorticities form a checkerboard 
matrix mg — ±^(1 -I- (—1)^1+^^). This phase is a spe- 
cific example of a staggered flux phase investigated in the 
High Temperature Superconductivity literature [25]. An 
explicit ground state configuration is shown in Fig. (2). 




FIG. 2. A ground state of the Fully Frustrated XY model 
in an explicit gauge. The arrows denote the rotor directions. 
The numbers (1 or 0) at the centers of the plaquettes de- 
note their topological charge. The solid lines denote anti- 
ferromagnetic bonds and the dashed lines denote ferromag- 
netic bonds. Any uniform global rotation will produce an- 
other ground state as will an interchange between sites on 
the two sublattices- a translation by one lattice constant. 



To distinguish between the two viable ground states 
and to highlight their Z2 symmetry and pinpoint the 
breaking of this symmetry a staggered magnetization 



M^iV(-ir- 



(27) 



is defined [24]. The current patterns within the two dif- 
ferent chiral ground states are rotated by tt/2 relative 
to each other. Domain walls separate the two checker- 
board ground states- these cannot be transformed into 
each other by a continuous rotation. The walls consist of 
links which separate two plaquettes with the same sign of 
the chirality. In some numerical studies the onset of XY 
symmetry "breaking" is extracted from an analysis of the 
helicity modulus. As mentioned earlier, the helicity mod- 
ulus exhibits a universal jump at the Kosterlitz-Thouless 
transition [26]. The helicity modulus may be measured 
by applying strains on the Josephson currents [20]. Let 
us next consider a 2 x 2 block with / = 1/2 everywhere. 
Such a small fragment is depicted in Fig. (3) where, in the 
chosen gauge, A25 = A^g = tt with all other Aij = 0. 
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FIG. 3. A 2 x2 block with / = 1/2: A25 = A58 = tt. 



Within the open boundary condition minimum of this 
small 2x2 block: (hoi ^ 6t4 
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"89 = 
'/'ge = (t>36 = (t>23 = '^°^^^('^) ^^^ ^i5 = (t>52 = '/'es = 
058 — cos^^(-^). This state lies in the topological sector 

(mi — 0,mii — I, mil I — ^Tfniv — 0). A similar state 
appears in the sector (1,0,1,0). The energy of these 



states is eo)i"o,i 



-1,0,1,0 



-4v5. The energy per bond 
in this state is lower than that of the global (entire lattice) 
ground state minimum. This is hardly surprising: open 
boundary conditions were placed on this small 2x2 block. 
At worst the energy of this block would have coincided 
with the global ground state energy. 

Let us denote the amount by which the open bound- 
ary condition state improves on the global ground state 
configuration (as placed on this small 2x2 block) by 



A_ = 4\/5 - 6\/2 = 0.458991. 



(28) 



The minimum of all (1,1,0,0) or (0,0,1,1) states is 



-1,1,0,0 



e™o"i 1 = ^8 (corresponding to 025 



f58 



with all other gauge invariant bond angles 0^^) = 0). 
The lowest energies of all other non checkerboard topo- 
logical sectors are even more elevated than those of the 



(1,1,0,0) or (0,0,1,1) sectors. Thus the energy of all 
non-checkerboard configurations are elevated by at least 

A+ = 6\/2 -8 = 0.485281... (29) 

relative to the global ground state configuration. 



V. Za SYMMETRY BREAKING IN THE FULLY 
FRUSTRATED MODEL 

All reported here and in the sections to follow is new. 
We described the ground states of the Fully Frustrated 
XY model on the square lattice. The topological charges 
behaved as they should- repelling each other and in the 
process generating a checkerboard pattern. We now 
prove that the uniformly frustrated XY model has a Z2 
phase transition: a claim that, till now, has only been ver- 
ified numerically and suggested by approximate schemes 
[24] . We show that the existence of a finite temperature 
phase transition is a consequence of finite energy domain 
walls between the two ground state domains. Very in- 
teresting earlier works examining topological excitations 
(including domain walls) are found in [27]. 



A. Energy Gaps on the Square Lattice 

The topological charge sectors' minima are separated 
from each other by finite energy gaps. For the 2x2 
block in (Fig. IV) the lowest energy configuration which 
does not belong to the "good" sectors \rni ~ l,mii = 
0,mjii = l,miv = 0} or {mi = 0,to// = l,?7i/// = 
0, mjv = 1} topological charge sectors is given by (f>i4^ — 

047 = '/'78 = </'89 = 096 = 063 = 032 = and 025 = 058 = 

TT and it lies in {mj = 0,m// = O^mjjj = l,mjv = !}• 
The energy of this state is E^q\ i ~ —8. Similarly if 
025 — 058 = e — TT (with e = 0+) a state of a similar 
energy lies in {mj = l,mii = l,?7i/// = 0,?Ti/y — 0}. 
^110 ~ ~^- Thus the infimum of the energy in all 
other "bad" topological charge sectors is e^"^ = —8. The 
{0ij = ±7r/4} states, although they arc the global min- 
ima configurations, are not the lowest energy configura- 
tion on this fragment of the lattice with open boundary 
conditions. Nonetheless, these reference states lie in the 
"good sectors" {mj ~ 0,mji = l,mjij = 0,mjv = 1} 
and the {mj = l,mii = 0,miii = l,miv = 0}. The 
energy of this state on a 2 x 2 block is e^ q -^ q = — 6\/2. 



In a system with periodic boundary conditions 



E = 



1 



J^c(J^c 



1 



E 

Blocks 



e\ 



(30) 



where e\ is the net energy of all bonds lying in an Lc x 
Lc block A. The sum is performed over all overlapping 
Lc X Lc blocks. The multiplicative factor in front stems 



from the fact that each bond will be included in the sum 
Lc{Lc + 1) times if one were to cover a square lattice, 
with periodic boundary conditions, with all overlapping 
Lc X Lc blocks. We foucus on the simplest case: Lc ~ 2. 
In Eq.(29), A+ = e^^l^ - e^^/^o, with e^/ the energy of 
a reference global ground state configuration (i.e. {(j)ij — 
±7r/4}) when evaluated on a finite Lc x Lc plaquette 
block. Within each small block A, for states of the "bad" 
incorrect topological charge registry: 



A+<6^[{m.,.J^i(l + (-l)-^+-^) 



or -(l + (-l)-^+-^+^)]- 



c'-e/ 



(31) 



where, as before, (xi,a;2) label the plaquette coordinates 
in the plane. Averaging the last set of inequalities over 
the entire lattice (A): 



E%m^,^,}]>E, 



A, 



Lc{Lc + 1 



-Nh- 



Lc{Lc + 1) 



N, 



good 



(32) 



where E^ is the global ground state energy (of the lat- 
tice with periodic boundary conditions), Nb denotes the 
number of blocks with incorrect {mxix2} registry, while 
Ngood denotes the number of 2 x 2 blocks with the correct 
checkerboard registry. Henceforth, we denote 



A = A, 



A_ >0. 



(33) 



Similarly, for the unfrustratcd (/ = 0) system, an |to| = 1 
plaquette {Lc = 1) with open b.c. is elevated by, at least, 
A = ^ relative to the m — ground state configuration. 
Here we employ the last relation to obtain a lower bound 
on the energy gap given the number of blocks (plaquettcs 
in the / = case) which do not have the correct topolog- 
ical charge registry (i.e. have nonzero vorticity m ^ 0). 
The reader can see how these bounds may be general- 
ized to higher order commensurate values of / where the 
block size Lc will be larger. These simple inequalities can 
be seen as a generalization of the Peierls bound [28] . 

We are now ready to apply the machinery of reflection 
positivity. Let us quickly outline the underlying logic: 
We assume that at infinity we have "good" (72 bound- 
ary conditions, by which we mean that one of the two 
types of topological charge checkerboard configurations 
will be found. We consequently bound the size of the re- 
gion that is occupied by the other "good" checkerboard 
ground state 171. In essence, we bound the size of the 
region separating the two checkerboard configurations; 
in this region the topological charge configurations are 
"bad" and costly in energy. The number iVf, of the blocks 
which have "bad" topological charge configurations is, 
at least, linear in the perimeter of the regions {domain 
walls) surrounding the gi droplets. The correspondence 
to the Peierls proof of a phase transition in an Ising ferro- 
magnet becomes transparent toward the end of the proof. 



We avoid Peierls' contour inversion [28], opting instead 
for a direct computation of the ratio of partition func- 
tions in order to find the probabilities of various contours. 
We prove, trivially, that at sufficiently low temperatures 
the energy bounds on the cost of "bad" regions (domain 
walls) translate into bounds on their probabilities. 



B. Reflection Positivity 

We employ Reflection Positivity [29] as applied to a 
system invariant under a reflection about an axis P pass- 
ing through lattice sites. The uninitiated reader is re- 
ferred to Appendix A for a brief exposure to this tech- 
nique. We partition the plane into two sides which flank 
the axis P (which itself lies along a line of bonds): P±. 
Let us define V± as the union of P± with P. Under 
reflection about P: 



us next reflect bk to cover the entire lattice A and gen- 
erate a vortex pattern Cfc. One may define a partition 
function (Zc^. ) that is the sum of all of the Boltzniann 
weights corresponding to these configurations. There are 
two correct (checkerboard) topological charge configura- 
tions within the Lc x Lc block; we shall denote these topo- 
logical charge configuration by gi and g2- For all faulty 
(non gi or 92) blocks bk reflections generate Nb > t'N 
faulty Lc X Lc blocks within the lattice with t' > 0. As 
before, let {A^} = vortex patterns in the Lc x Lc blocks 
centered about the sites i, and let {171,52} denote the 
two correct (checkerboard) vortex configuration within 
an Lc X Lc block. Note that the set {Xi} consists of all 
maximally overlapping Lc x Lc blocks. 



C. The Kaleidoscope Patterns: Energy Gaps on The 
Reflected Blocks 



RpF{S{xi), ..., S(xn)) = F{S{RpXi), ..., 5'(i?pf„)) (34) 
with Xi planar locations. Thus, under reflection 



^p^ij 



RpOi 



if (ij) is vertical: Rp4>ij 
if (ij) is horizontal: Rp4>ij = 



ORi 

hij. 



(35) 



For both horizontal and vertical bonds, cos 0ij is invari- 
ant under Rp. Defining 

B = -[ Y^ -- ^ ]cos(/)ij, 
{ij)ev± {ij)eP 



(36) 



the Haniiltonian H = B + RpB is manifestly invariant 
under Rp (yet the reflected plaquettes correspond to the 
same frustration < / < 1/2 only when f ~ 0, 1/2). As 
the sense of the circulation in the sum X]n_R '/'y i^ fixed 
over the entire lattice, under Rp-. {mg — f) -^ —{mg — f). 
For / = ^, the vorticities transform as 



-1,0,1,2^2,1,0,-1 



(37) 



under reflections. The correct (ground state pattern) set 
of topological charges {mg\ at the various planar loca- 
tions X is invariant under Rp. This is to be expected 
as the energies are unchanged by Rp- there is no energy 
increase along the reflection axis P. When / = 0, the 
vorticities transform as 



-1,0,1^ 1,0,-1 



(38) 



under Rp. For < / < 1/2, the quantity 27r[— (ttt,- /)] is 
not an allowed value for the circulation sum. Let us de- 
flne bk as the Lc x Lc topological charge pattern about the 
central plaquette at Xk, and introduce a function uji,^{x) 
which is equal to one if bk is found within the 2x2 block 
centered about x, and is equal to zero otherwise. Let 



By repeatedly reflecting any Lc x Lc block A to cover 
the entire lattice we generate a periodic pattern of repete- 
tive 2Lc X 2Lc super-blocks. The planes P may be 
either vertical or horizontal and for each orientation 

Rl hor^zontal = ^p vevHcai = "^- After two reflections by 
any two parallel planes the original pattern must emerge. 
Taking note of all possible "bad" (non-checkerboard) 
Lc X Lc blocks one sees that in all cases the global pattern 
Ck generated by consecutive reflections contains more 
bad Lc X Lc overlapping blocks than those of the checker- 
board type: in Eq.(32) iVf, > Ngood- Consequently, the 
energies of the all repetitive patterns satisfy 



Eck > Eq 



>E^ 



A 



Lc(Lc + l 

A 



-m 



Lc{Lc + l) 



t'N, 



(39) 



with Lc = 2, and, as denoted earlier, with the number of 
bad blocks in each repetitive pattern Ni, > t'N. 

D. Chessboard Estimates 

Let us next deflne an angular volume Q\ within a LcX 
Lc block such that sup^gj^A 1^' ^ Qaroun 



i (where 9f 



round 



I = S for all sites 
denotes the angular configuration of a 
global ground state minimum on the entire lattice- the 
pinned down reference angular configuration with broken 
global Z2 and U{1) symmetries) and, correspondingly. 



max ^ ^ supJea} = Cre/ + eo- 

LaxLahlock partitions a 



(40) 



The size of il\ is \il\\ = (26)^. The volume f^A is strictly 
contained in a single checkerboard sector. On the fringes 
(dftx) of this volume: 



e\dn, < e^'^f + eo 



(41) 



with e^'^f the energy of a global ground state configura- 
tion when evaluated on a finite Lc x Lc plaquette block, 
and eo < A. Next, we define a global angular volume fi 
such that supggfj \6i — Qsroun I _ J |.^^ ^j^ lattice sites i. 
The energy of any global lattice configuration within this 
volume has an energy which is bounded from above by 



Edn < Eq 



Neo 



Lc{Lc + l) 



(42) 



The definition of 51 and Q\ is made possible by the 
continuity in {di} of the Hamiltonian in Eq.(l). Now 
[Zc^^/Ztot] may be easily bounded at low temperatures 



Ztot Zn 






/3t'N 



A- 



eo 



LciLc + 1) 



(43) 



where (27r)^ is the net angular volume of the entire sys- 
tem and {25)^ is the volume of il. Here we noted the 
angular volume of Ck is bounded by {2tt)^ . [30] 

Let us define a "logical" test function on a block \i 
centered about site i by F{Xi = (71,32) = 1 and F{Xi ^ 
91.2) = (the reader will note that is merely a variant of 
Wfij, (x) : F = LOg-^ + ujg^ ) . Let us furthermore define 



vix) ^ n ^(^^ 



(44) 



Se^i 



The region where V{x) vanishes defines a domain D. In 
the spirit of providing very generous upper bounds, let i 
denote the outer perimeter of D, and let us pick merely 



f ^max{Int[^],l}, 



(45) 



non-overlapping faulty bk (where Int[ ] denotes the inte- 
ger part). The probability of a domain D, 



Prob{D) < max(a;bj(xi)...cj6^(x£)) 



(46) 



fc=i 



with 



t = [t'(A - eo) - r^ Htt/S)]/{L,{L, + 1)). (47) 

In the last inequality, the blocks {A^} are chosen to be 
non-overlapping: no spins are shared by two blocks (in- 
cluding the boundaries of the blocks). The second in- 
equality in Eq.(46) followed from reflection positivity 
(The reader is referred to equation (70) in Appendix A 
for details.) Let G* {i = 1,2) denote the support of 
the correct vortex configuration (gi) in a region of any 
size. Let us assume that at infinity, we have G^ bound- 
ary conditions. The probability of having a domain D 
incurred by Nqi units (in a partitioning of the lattice 



to non overlapping blocks) of the opposite checkerboard 
state is bounded by Eq.(46). The bounds on the area 
covered by the opposite checkerboard state may proceed 
as in the usual Peierls argument (see Appendix B and 
Eq.(74) for details), to yield 

.N, 



^-(T)< E ^^3^ 



=4Lc 



VI. 



X ex.p[~l3Wt max{Int[^], 1}]. 



(48) 



The right hand side can made as small as desired for 
/3 > Pcir). Stronger bounds hold for all other "bad" 
topological charge configurations. Thus, at sufficiently 
low temperatures, a spontaneous symmetry breaking is 
sparked by applying boundary conditions at infinity. 

Analogously, for f — 0, the fraction of |?7i| ~ 1 vortices 
(which, for Lc = 1, are elevated by at least A ~ 7/2 rel- 
ative to the lowest lying m ~ configuration)is bounded 
from above by arbitrary, exponentially small, numbers at 
sufficiently low temperatures: "pair annihilation" . 

Here we avoid the contour flipping algorithm of the 
standard Peierls argument by use of Reflection Positivity 
to bound Prob{D) in Eq.(46). We kept Lc general (and 
have not explicitly set Lc = 2) in order to suggest other 
plausible extensions (not hinging on Reflection Positiv- 
ity) to other values of / for which Lc will be larger). To 
emphasize- in this proof a necessary (but not sufficient 
requisite) is to establish a gap between the minimum en- 
ergy per site of "bad" topological configurations on the 
Lc X Lc block with open boundary conditions and thermo- 
dynamic limit ground state tiling the entire plane belong- 
ing to the "good" topological sector (which we dubbed 
the "reference" state). In general, the open boundary 
condition minimum on A will cannot tile the entire plane 
and will produce a slightly lower ground state energy per 
site than the ground state of the entire lattice. 



VI. AN ALTERNATIVE PROOF 

The reader might be a bit dismayed that our proof 
seems to rest on the relatively fortuitous event that for 
the small Lc — 2 the gap A = A-|_ — A„ > 0. Had this not 
been the case we might have been required to compute 
energies on much larger Lc x Lc blocks. Moreover, in the 
previous proof, in subsection(VC), an explicit check that 
the number of bad blocks N^ is greater than the number 
of good blocks was necessary for the configurations gener- 
ated by consecutive reflections of {bk} to cover the entire 
lattice. Here we show that this is not the case when at- 
tacking the problem along a slightly disparate path. In 
the spirit of Eqs. (36,57) let us define 



EA' 



E< 

{ij)eR 



- E 

2 ^ 



cos (f>i 



(49) 



for each s x s block A' (with an arbitrary integer s > 1) 
such that 



A' 



(50) 



in a partition of the lattice with "marginally" non over- 
lapping blocks A' (nearest neighbor blocks share only a 
common boundary dX'). The form Eq.(49) is particularly 
suited for Reflection Positivity treatments. Note that the 
open boundary condition ground state {6i} of e\i on the 
block A' (the previous "reference state" ) may be repeat- 
edly reflected (as / = 1/2 is congruent to (— /)) along the 
boundaries 9A' to tile the entire plane. Thus the open 
boundary condition minimum on the small block is, in 
this case, a portion of the global ground state. 

Now, a flnite energy gap A > (necessitated by our 
deflnition of "bad" and "good" topological sectors) sepa- 
rates the various "bad" topological charge conflgurations 
(which do not contain a ground state of e^i on the fl- 
nite s X s slab with open boundary conditions) and the 
ground state value of e^' . The reader might recognize A 
as none other than A+ for the particular case s = Lc as 
deflned for the energy in Eq.(49). In computing the prob- 
ability of a domain D we note that net domain size A^^ 
depends on how we partition the plane into marginally 
non overlapping blocks A'. If, in a partitioning of a lattice 
into maximally overlapping blocks {A}, we find A^^ faulty 
blocks of the incorrect topological charge registry then in 
the optimal partitioning of the lattice into marginally non 
overlapping blocks {A'} we will find, at least, 



iV^ = Int 



Nh 



s^-1 



(51) 



faulty blocks {AJ,}. Consequently, in Eq.(47), we may 
substitute 



A 



A 

l2 



(52) 



and repeat the previous arguments of section (V). 

Although, in the above, we focused attention on s x s 
square blocks, we could have chosen A' to be a domino: 
the minimal (2 x 1) plaquette configuration. There are 
two checkerboard 171^2 states (the sectors {mi, 1112) — 
(1,0) or (0,1)) within the domino and all of our argu- 
ments may be reproduced with the factor s^ of Eqs.(51- 
52) replaced by the area of the domino (two plaquettes). 



VII. THE STANDARD XY MODEL- A PROOF OF 
VORTEX CONFINEMENT 

We may similarly bound the number of radius r vor- 
tices (r 3> 1), which were defined previously, within the 
unfrustrated system. When r ^ 1 the energy penalty is 
no longer a mere 7/2 J as it was for a vortex of size r = 1 



(which was immediately adjacent to an antivortex), but 
rather a forbidding logarithmic function of r. Employ- 
ing the Reflection Positivity inequality Eq.(65) with the 
energy penalty associated with the large radius vortex of 
Eq.(25), the probability of a vortex of size r 



Prob{r - vortex) < {n/rf^^' 



(53) 



where, we once again employed the fact that the volume 
of the state generated by repeated reflections of the ra- 
dius r vortex to cover the entire plane has an angular 
volume which is trivially bounded by {2tt)^ . This prob- 
ability can be made as small as desired at sufficiently low 
temperatures (large /?). The probability of a radius r 
vortex drops algebraically with r: There are no "real" 
macroscopic vortices at sufficiently low temperatures. 



VIII. KALEIDOSCOPE FIELD PATTERNS 

Although in many instances (e.g. uniform non com- 
mensurate /) the determination of the exact ground state 
is nontrivial there are several exceptions to the rule. 

If the frustrations {/a'} form a kaleidoscope pattern- 
i.e. if they may be generated, starting from an arbi- 
trary sequence of frustrations in a small region, by con- 
secutive reflections in different planes to cover the en- 
tire lattice, then a ground state may also be formed by 
the same sequence of reflections. Let us now "elaborate 
" and prove this statement. Let X'l be a block span- 
ning one or more plaquettes on which the frustrations 
at the different plaquettes are arbitrary. If the block is 
of flnite size then the determination of the exact open 
boundary condition minimum of the energy functional 
ey deflned in Eq.(49) will yield an angular conflgura- 
tion(s) {9i} on X'l. Now let us invert the frustrations 
and set fy — —fy for the block A2 which is the mir- 
ror image of A'^ about a plane Pi which passes through 
one of the edges of A'^ : X'2 = Rp^X'^. We now note 
that the open boundary condition minimum of ey on the 
fused super-block A'^ U A2 must satisfy min{eA' ua' } = 
minJEAi + ^x^} ^ [minjeAi} + minjeAa}], Y^t just such 
a conflguration saturating the lower bound with an en- 
ergy [minjeAil + rnin{eA2}] may be constructed. To do 
so we may simply leave the original angular conflgura- 
tion {6i\ for all sites i within the block A'^ untouched 
and set 6i = 6jip i for all sites i in the reflected block 
A2. We may now choose a different plane P2 and keep 
repeating the process until A'^ U Aj U A3 U . . . spans the 
entire lattice. Perhaps the simplest realizations of such 
ground states are the staggered flux phases (of which the 
/ = 1/2 ground state hitherto considered is a special ex- 
ample). As other realizations we may consider systems 
in which the frustration 



is = I sm( ) sm( ) 

n n 



(54) 



10 



for all plaquettes x within the dual lattice with an arbi- 
trary f and a finite integer n. Here X'l is the elementary 
n X n block. Any reflection of the frustrations fg^y, in- 
verts their sign on the sites RpX which lie in A^_|_]^. Many 
systems, such as that of Eq.(54), are reflection positive 
and a proof of a finite temperature ordering similar to 
that for / = 1/2 follows. 



IX. GENERALIZED Zk SYMMETRY AND 
MAGNETIC GROUPS 

For a uniform frustration / = p/q with relatively prime 
p,q > 3, the topological charge ground state consists of 
repetitive units extending ki,k2 units along the a;i,a;2 
axis respectively. By translating the origin in the cover- 
ing of the lattice by non overlapping blocks, we may gen- 
erate k — kik2 ground state configurations within the 
repetitive blocks. To describe this viable translational 
Zfej (g) Zk2 symmetry, we may extend the definition of the 
Z2 order parameter previously defined: 



^k = im]^s(^) cxp[27ng/k]l 



(55) 



where the sum is now over non overlapping blocks A and 
possible ground state topological charge sectors g. As be- 
fore, i-Ug(X) = 1 if the pattern bg appears in A and is zero 
otherwise. We may similarly define the two block corre- 
lator {uji(X)LUj(y)) for two blocks A and A' appearing in a 
partitioning of the lattice into non overlapping blocks. In 
the limit q ^ 00 (irrational /) the discrete translational 
symmetry transforms into a continuous symmetry 



Z, 



ki 



0(2), 



(56) 



and the translational symmetry group becomes the sym- 
metry group of the two-torus. 

The Z2 and 0(2) symmetries present in the / = 1/2 
case become two 0{2) symmetries in the limit of large q 
(continuous /) . Physically, all possible translations of a 
given ground state configuration amount to a rotation of 
the spins by an all possible angles. One might speculate 
that perhaps for irrational /, there is no finite temper- 
ature vortex ordering (no explicit Zfe. -^ 0{2) symme- 
try breaking) but rather a finite temperature Kosterlitz- 
Thouless like transition (in which the vortex ordering "al- 
most" breaks an 0(2) translational symmetry but not 
quite) [34]. Such a viable link between the degeneracy 
amongst states generated by displacements and those cre- 
ated by global rotations is also found elsewhere. In the 
Quantum Hall problem, the degenerate states may be 
labeled in terms of their angular momentum quantum 
numbers or, alternatively, by their magnetic translation 
operator eigenvalues. On a lattice (the Hofstader prob- 
lem) the correct degeneracy is easily read off in terms of 
the magnetic translation operators; the flux piercing the 



fundamental cell is an integer multiple of the elementary 
fluxon. An irrational fiux / might lead to an effective 
loss of commensurability effects (the size of the elemen- 
tary cell becomes of the order of the system size) and 
the equivalence between those degeneracies spawned by 
translations and those created by rotations may be re- 
stored. The sole (albeit important) distinction between 
the Quantum Hall problem and the frustrated XY model 
is that in the latter the magnitude of the wavefunction is 
pinned, |V'(2i')| = 1, whereas in the Quantum Hall prob- 
lem both phase and amplitude vary. 

In addition to the trivial Zk^ groups generated by dis- 
placements, one finds in certain instances a greatly en- 
hanced symmetry group. The ground states of the frus- 
trations / = 1/3 and / = 1/4 are a case in point [37]. 
Here a multitude of "zero-energy domain walls" appear. 
Phrased in our language this implies that nontrivial cov- 
erings of the lattice are possible such that within all over- 
lapping blocks A we find the reference global ground state 
minimum configuration. At low temperatures the infinite 
degeneracy is lifted by the free energy of spin waves. 



X. ISOTROPIC TO NEMATIC TRANSITIONS 

Just as we were able to bound, at low temperatures, 
the number of free vortices within the unfrustrated XY 
model we may also prove that free disclinations in liquid 
crystal system must condense at sufficiently low temper- 
ature; such a confinement of disclinations occurs in the 
transition between the isotropic and nematic (or "topo- 
logically" ordered [31]) phases. Similar arguments apply 
for other systems having other order parameter symme- 
tries. Within the (unfrustrated) XY model, the 0(2) 
vortices are elevated by a finite amount relative to the 
ground state. Similarly within lattice models of liquid 
crystals the RP^ topological defects (disclinations) of the 
headless director fields W must incur an energy cost rela- 
tive to the defect free ground state. This fused with a Re- 
flection Positivity symmetry which is present in many of 
these models (for some lattice models sec [31,32]) allows 
one to set bounds, via chessboard estimates, on the fre- 
quency of disclinations at sufficiently low temperatures. 
Analogously, one may examine an analogue of the Fully 
Frustrated XY model: a nematic on a lattice having half 
an RP^ fluxon thread each of its plaquettes. Similar ar- 
guments may be used to bound the frequency and size of 
topological defects associated with different systems. 

The chessboard estimates might find a natural arena 
in many of the lattice gauge theories. For instance, lower, 
area law, bounds on the Wilson loops in Z2 lattice gauge 
models may be immediately derived- no high tempera- 
ture expansions are necessary to prove this point. These 
bounds follow effortlessly from Reflection Positivity in- 
equalities, Eqs. (64,65,70): The Wilson loop is the quan- 
tity appearing on the right hand side of these inequalities 
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and the identity a"^ = 1 {A = Yl^j a) may be employed 
for the repeated (reflected) elements of Z2 appearing on 
the right hand side to prove a lower, area law, bound. 



XI. CONCLUSIONS 

We presented a theoretical study of the XY model with 
specific emphasis on the fully frustrated case. The main 
result is a proof of the existence of a Z2 symmetry break- 
ing transition within the fully frustrated XY model con- 
sistent with previous numerical results available in the 
literature. We further discussed extensions of the meth- 
ods employed to other arenas and models (confinement of 
disclinations in the isotropic to nematic transition, topo- 
logical charge orders in various XY models and bounds 
on vortices within the unfrustrated XY model). 
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XII. APPENDIX A: WHAT IS REFLECTION 
POSITIVITY? 

Reflection Positivity has long been a weapon reserved 
for the cabinets and drawers of MPs (Mathematical 
Physicists). Our simple minded nuts and bolts approach 
presented below is geared toward our very specific limited 
set of problems and is quite superficial. It does the field 
no justice. For a more detailed and rigorous exposition 
the reader is referred to [29]. 

Let us partition the plane into two sides which flank 
an orthogonal bisecting plane P (which itself lies along a 
line of bonds): P±. Define V± as the union oi P± with P. 
In a Reflection Positive system, there exists an operator 
B such that the Hamiltonian may be expressed as 



H = B + RpB. 



(57) 



As noted earlier in the text, the Reflection operator Rp 
about the plane P is defined by 

RpF{S{xi), ..., S{xn)) = F{S{RpXi), ..., S{RpXn)) (58) 

where Xi denotes the planar location of site i. [33] 
The basic inequality: The thermal average 



{ORpO) > 



(59) 



for all operators {O}. (In other words, such thermal 
averages behave like inner products in a "positive" sense.) 
Proof: Let us write the thermal average longhand 



{ORpO)^Z-^ E E E ^^p^ 

{s^}eP+ {Si}<£P- {s^}eP 

exp[-f3{B + RpB)], (60) 

where the partition function Z > 0. Thus the sign of 
the thermal average is sign of the multiple sum which, 
by symmetry, may be folded back onto one half plane 
to become a sum of squares when {O} is a local opera- 
tor defined on P+ or on 7-"-. Let F = Oexp[—PB] be 
an operator on V+. Let us explicitly segregate the spin 
arguments in V+ which correspond to spins lying on P 
(which wc shall denote by {S{xj)}) and those that lie in 
P+ (labeled by {Six^)}). Note that (Eq.(58)) RpF on 
V- is identical to F on V+. Thus Eq.(60) reads 

{s,}eP {s,}eP+ 
F{{S{x, e P+)}, {Six;, e P)})] ' > 0, (61) 

which proves the prinicipal Eq.(59). Inserting O = A — 
B{BRpA)/{BRpB) in Eq.(59), a Schwarz inequality 

{ARpBf < (ARpA) X (BRpB) (62) 

results. Setting O = A - (A) in Eq.(59) (or B = 1) 

(A) 2 < (ARpA). (63) 

We remark that this general inequality immediately im- 
plies that opposite topological charges always bind [35]- 
a similar yet weaker version of the relations derived in 
[36] . We may repeat the last step n times in a row to get 



{Ay < (A{Rp,A{Rp,A{...Rp^A)...))). 



(64) 



If we set the size of the system to be A^ = 2" then the 
last inequality would imply that 



{local A) < (Reflected A) 



(W/N) 



(65) 



where local A denotes the local function A and 
Reflected A denotes the function A generated by consec- 
utive reflections of A in different planes, until it covers 
the entire lattice. We have more cautiously inserted a 
constant W instead of one. The functions A that were of 
interest to us were local test functions that were equal to 
one if a specific local configuration was found and were 
zero otherwise. In this manner the probability of finding 
certain necessarily "bad" conflgurations, lying on the in- 
terface between the two ground state domains, may be 
bounded. Here the thermal average (A) is the probability 
that such a "bad" configuration would be found. 



(Reflected A) 



Zpeflected A 



(66) 



where Zpefiected A is the sum of all Boltzmann weights 
corresponding the Refiected A configuration which spans 
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the entire lattice, and Z is, as usual, the sum of all Boltz- 
mann weights- the partition function. Eqs. (65) and (66) 
bound the probabilities of certain local configurations. 

We will now fortify and extend these relations to reach 
the more sophisticated inequality (Eq.(70)) employed in 
our proof. Suppose that A^ denotes the observation of a 
certain local configuration about point x\ while A^ de- 
notes the test function for another configuration being 
seen at x^- By the Schwarz inequality (Eq.(62)) 

(AiA2)2 = (Aii?p32)2 < (^ii?pAi) X {A^R^A^), (67) 

where A = RpA. We now repeatedly apply Eq.(63) (as 
in Eq.(64)) to the right hand side of Eq.(67) to obtain 



{A1A2) < (Reflected Ai 



AW/N) 



X {Reflected A2)'^^/'^\ 



(68) 



which may be followed by an insertion of the identity 
in Eq.(66). For the correlation function of three local 
configuration test functions we may note that 



{AiA2A3f = {{AiA2)RpA3) > 
{{AiA3)RpA2) X {{A2A3)RpAi) 



(69) 



i.e. a product of terms of the form (ARpB). We may 
consequently apply the Schwarz inequality, as was done 
in Eq.(67), and then invoke Eq.(66) to get a simple ex- 
tension of Eq.(68). Thus, as the reader can see, for the 
i?-th order correlator we have 






eflected Ai 



z 



W/N 



(70) 



Eqs. (70) are a central result. With these relations (the 
chessboard estimates), local correlations (appearing on 
the left-hand side) may be efhciently bounded in terms 
of global events (which form the right-hand side). 



XIII. APPENDIX B: CONTOUR ARGUMENTS 

In a two dimensional nearest neighbor Ising system 
the probability of having a certain specific contour of 
length £ separating an island of "-" ( "Nq ) spins in a sea 
of "-I-" spins (or vice versa) is bounded by exp[—2pj£] 
when the exchange constant units are restored. Here we 
repeat a variant of the standard Peierls argument for the 
benefit of our readers. Let us assume that there are "-I-" 
conditions at infinity. We may evaluate the probability 
that the lattice site af = is occupied by a "— " spin. To 
have a minus spin at the origin there must be at least 
one contour F encircling the site 0. 



Prob{S{x = 0) = -1) < Yl Prob{r) 
r enc. 

The probability of a given contour F of length i: 



(71) 



ProbiD = f < -^i— 
"^ ^r flipped 



= exp[-2/3J£], (72) 



where Zp fjippcj is the sum of Boltzmann weights over 
all spin configurations in which the contour F is flipped: 
All of the spins inside F are flipped with all other spins 
in the system are untouched. The energy of each in- 
dividual configuration summed upon in Z-p is exactly 
2J£ higher than that configuration in which F is flipped- 
whence the second equality in Eq.(72) follows. A mo- 
ment's reflection reveals that Zp flJQpjgj only contains 
a subset of all the terms which are present in Z and 
consequently Zp flipped < Z, implying the inequality in 
Eq.(72). Eq.(71) implies that 

Prob{S{x = 0) = -l) < Y^ n(^)exp[-2/3J^]. (73) 

i=4,6,... 

We now bound the number, n{£), of contours of speci- 
fied length £. This number is bounded from above by the 
number of unrestricted random walks of length £. If given 
a unique starting point, the number of non-backtracking 
random walks of length £ on the square lattice is bounded 
by 4 X 3^^^. The number of possible starting points for 
a random walk enclosing the origin (x = 0) is strongly 
bounded by the area enclosed by the contour. For a given 
contour length £ this maximal area is (£/4)^, correspond- 
ing to a square of side (^/4). Bounding this maximal area 
by (3£^/4) provides a very generous over estimate. 
Fusing all of the bounds together: 



ProbiSix^O) = -1) = 



.Nc 



N 



< 



Y ^^3^exp[-2/3J^]. 



(74) 



£=4.6, 



The last sum can be made smaller than 1/2 (indicating 
symmetry breaking) for sufficiently large /? (low temper- 
atures). The proofs presented for the XY models (sub- 
section(VB)) rest on the same logic with the twist that 
the contour inversion probability bounds (Eq.(72)) are re- 
placed by the Reflection Positivity inequality of Eq. (70) . 
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A vortex may be described by a source term, 

H + AiV X V X _ff = qaz5'^{x) 
with go the flux. In cylindrical coordinates 



H, 



-v^-^^-^^-^^^'^^ 
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The solution with the required singular behavior near the 
origin is 



H,= 



10 y ( r s 



(82) 



with Ko a Bessel function. If there are two vortices then 
on the right hand side of the London equation we should 
insert two source terms: z[qi5^{x — xi) + q25^{x — X2]. 
The corresponding solution is a simple super-position. 



H, 



2-k\ 



Vbi^o(^ 



Xl\ 



Xi 



+ q2KQ{ 



■ X2 



Xi 



A^ Al al 

Inserting this in Eq.(78), the interaction energy 

— ==2(^^^^^)7^o(^^^). 
Asymptotically, 

Ko{x) - ln(i) (a; < 1) 



)]. (83) 
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Ko{x) ~ W — exp[-2;] {x > 1) 
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47rne2 



(V X Hfd^x (77) 



where we invoked Ampere's law VxH = —j. The energy 
(kinetic -|- field) of a vortex line of length L 



Stt 



d^x [H'^{x)+Xl{V xHf] 
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